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Abstract 

In this paper we show that the so-called array Frechet problem in 
Probability/Statistics is (max, -|-)-linear. The upper bound of Frechet 
is obtained using simple arguments from residuation theory and lattice 
distributivity. The lower bound is obtained as a loop invariant of a 
greedy algorithm. The algorithm is based on the max-plus linearity of 
the Frechet problem and the Monge property of bivariate distribution. 

Keywords: Max-plus algebra, Frechet bounds. 

1 Introduction 

As a preliminary remark the author would like this paper to be a modest 
tribute to the work of Maurice Frechet in Statistics. The work has been 
started at the occasion of his 130th birthday and the 100th anniversary of his 
stay in Nantes (the town the author is living in) as professor in Mathematics. 

In this paper it is shown that the tropical or max-plus semiring Kmax (i-e. 
the set R of real numbers with max as addition and -|- as multiplication, see 
the precise definition in subsection 12. 2p is the underlying algebraic structure 
which is well suited to the Frechet contingency (or correlation) array problem 
[3]. In other words the Frechet problem is a tropical problem which thus 
has its place in the new trends of idempotent mathematics founded by V. P. 
Maslov and its collaborators in the 1980s (see e.g.[7| and references therein). 

Prom this main result the Frechet upper bound is derived by residuation 
and the distributivity property of Mmax as a lattice. The Frechet lower bound 
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is obtained as a loop invariant of a greedy algorithm. This algorithm is based 
on the tropical nature of the Frechet problem and the Monge property of 
bivariate distribution. 



1.1 Organization of the paper 

The paper is written to be sef-contained. Thus, in Section [2] we introduce 
main notations used in the paper, we restate the Frechet array problem and 
its bounds. We define the tropical semiring Mmax and recall basic results 
on residuation theory and lattices. In Section [3] we prove the main result of 
the paper that is the Frechet array problem is max-plus linear in the space 
of cumulative distribution functions (see Theorem 13. From this result in 
Section [J] we derive the Frechet bounds using new approaches. The upper 
bound is derived from residuation theory and the lattice distributivity prop- 
erty of the max-plus semiring Mmax (see Corollarv 14. ip . The lower bound is 
obtained as the loop invariant of a greedy algorithm (see Proposition 14. 1 p . 
In Section [5] we conclude this work. 

2 Preliminaries 

In this Section we recall basic results concerning Frechet array problem and 
the max-plus semi-ring Mmax- 

2.1 The Frechet contingency array problem and its solution 

This problem is described in e.g. Let n be an integer > 1. The set 
Matnm(M+) denotes the set olnxm matrices which entries are nonnegative 

V 

real numbers. We define the partial order < on Mat„„(R_|_) as follows: 



A = [ai^j] ^ B = [bij] 4> Vi, j, ^ ^ ai^k < X] X] 



(1) 



1=1 k=l 1=1 k=l 



Introducing the fundamental n x n matrix: 




(2) 



V 

where l{i<j} = 1 if ^ < j and otherwise, the partial order ^ can be 
rewritten as follows: 



A^B ^ DAD^ < DBD^ (entrywise) (3) 
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where () denotes the transpose operator. 

Let p,q G Mat„,i(]R+) such that J2^^iPi = = X]j=i 1j- Without loss of 
generahty we can assume: 

0- = 1. 

The problem of Frechet is then to find (if exist) the maximum and the 

V 

minimum element w.r.t the partial order ^ of the subset of Mat„„(IR+): 

n{p, q) = {F G Matnn(K+)|i^ satisfies (Fl) and (F2)} (4) 

with: 

(Fl). Fl=p, 
(F2). l^F = qT. 

Where 1 denotes the n-dimensional vector which coordinates are I's. 

Frechet proved that there exist a maximum element, V'd ^(Pi o) "= 

^max5 and a minimum element, ^{{p^q) Fmin such that: 

{DF^^^D'^^, = mm{{Dp)^, {q^D^),) "= (F^ax)*j, (5a) 

and 

(DF^inD^),,,- = max(0, (Dp), + {q^D^)j - a) =' (F^in)ij (5b) 
for all i, j = 1, . . . , n. 

2.2 The max-plus semiring Mmax 

Let M be the field of real numbers. The max-plus semiring Mmax is the set 
M U {— cxd} equipped with the addition © : (a, 6) i-^ a 6 =^ max(a, h) and 

def 

the multiplication : (a, 6) i-^ a 6 = a + 6. The neutral element for is 
O := — oo and the neutral element for is 1 := 0. The max-plus semiring is 
said to be an idempotent semiring because the addition is idempotent, i.e. 
a a = a, Va. 

The idempotent semiring Mj^^x idempotent semirings isomorphic to 
Mmax has many applications in discrete mathematics, algebraic geometry, 
computer science, computer languages, linguistic problems, optimization 
theory, discrete event systems, fuzzy theory (see e.g. [1], [3], [6], [8]). 
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2.3 Order properties of Mmax and residuation 

Let us consider the max-plus semiring I^max cilrea^dy defined, in the intro- 

def 

duction. The binary relation TZ defined by: alZh ^ a ® b = b coincides 
with the standard order < on M. We denote M^ax the semiring completed 
by adjoining to Mmax a T := +00 element which satisfies a © T = T, Va, 
a T = T, Va 7^ o, and o T = o. This is mainly motivated by the fact 
that some of the further results can be stated in a simpler wa.y in Mj^-^ax- 
The completed max-plus semiring Mmax is a complete sup-semilattice, i.e. 

— Hpf — Hpf — 

C M ©yl = © X exists in M = M U {-cx), +00}. This implies that Mmax 

is also a complete inf-semilattice because /\A = ©{x G M|Va ^ A, x < a}. 
Thus, Mmax is a complete lattice. Finally, let us mention that Mmax is a 
distributive lattice, i.e.: 

o© (6 Ac) = (o A6) © (a Ac). (6) 

As we will see in the sequel this property will be of particular importance. 

— def — 

Let us define left and right division in Mmax by: h/a = ©{x G M|x a < 6} 

def — 

and a\b = ©{x G M|a0x < 6}. Left and right division coincide with 
the usual subtraction to which we add the following properties: o\a = T, 
T\a = o if a 7^ T, T otherwise (similar formulae for /). 

We extend operations and binary relations from scalars to matrices as 
follows. If A = [ajj], B = then: A< B (entrywise) 44> Vi,j, ajj < ftjj, 
A(BB = [ojj © 6jj], A/\B = [ojj A and AqB denotes the matrix such 
that its entry (i, j)i © is: {A B)i^j = (BkO-i^k © bkj- We obviously 

have: (^4 B)^ = B^ . We also extend the divisions to (possibly 
rectangular) matrices with suitable dimension: 

Hpf 

{A\B\, = {®{X\A QX< B})ij = Akak,i bkj, (7a) 

Hpf 

{D/C\j = {®{X\X QC< D}\j = Aid,,i/cj,i. (7b) 

The © in the formulae (j7ap and ()7bp corresponds to the supremum w.r.t 
entrywise order between matrices. The application Y A\Y (resp. Y 
Y/C) is called the residuated mapping of the application X Aq X (resp. 
X ^XQC). 

For more details on max-plus algebra and residuation theory we refer 
the reader to e.g [U Chp. 4 and references therein]. 
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3 Main Result 



We begin this Section by the following fundamental lemma. 

Lemma 3.1 Let {u^)Y=i m elements o/Mat„i(M+). Define = X]i=i '^^ ■ 
Then, 

[n^-.-n"] 1 = [s^---s"]0l (8) 

where 1 denotes the n- dimensional vector which all components are 1 's. 

Proof. Let us remark that: [u^ ■ ■ ■ u"] 1^ = s". Because the vectors have 
all their coordinates nonnegative: < ■■■ < s" (componentwise), which 
is equivalent to: s" = © • • • © s"". Now, we just have to remark that: 
© ■ ■ ■ © s" = [s^ ■ ■ ■ s^]Q 1 which ends the proof of the result. □ 

Remark 3.1 In the previous Lemma define U = [u^ ■ ■ ■ u"] and S = [s^ ■ ■ ■ 
We remark that S = UD^ , recalling that D is the matrix defined by 
Thus, relation (0) can be rewritten: 

Ul = {UD^)Q1. (9) 

Theorem 3.1 (Main Result) Let us consider a matrix F € Matnn(M+)- 
Let p and q be two elements of Matni{^+)- Then, 



j,^^/ {DFD^)Q1 =Dp 



Proof. Fenip,q)^{ f^"^ 



T 



IF =q' 

Because matrix D is invertible so is D^ and: 

Fl =p ( DFl = Dp (eq 1) 

l^F =qT ^\ iTpD^ = q^D^ (eq 2). 

For (eq 1). We apply Lemma [3.11 with := {DF)^j and = {DFD^) j be 
the jth column vectors of matrices DF and DFD^, respectively. We obtain 
the following equality: Dp = DFl = {DFD'^) 1. 

For (eq 2). We apply LemmaOwith := (-DF^).,j and s^ := {DF'^D'^),j. 
We have: DF'^1 = (g^L>^)^ = (DF'^D^) 01. By definition of ()^ we have: 
q'^D'^ = 1^ {DFD'^) which ends the proof of the Theorem. □ 
This result can be reformulated as follows. Let us define the following 
sets: 

Hpf 

Dist„,i = D Mat„i(M+) = {Dx, x € Mat„i(M+)} (10) 
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and 

DisW = D Mat„„,(R+) = {DXD' , X G Mat„„(M+)} (11) 
and for all P,Q Dist„i: 

H(P,Q) = (F £ Dist„„|F0l = Pand 1^0F = g^}. (12) 
Then, Theorem O states that Vp, g G Mat„i(M+) and VF G Mat„„(R+): 

F G n{p, q) ^ DFD^ G m{Dp, Dq). (13) 

Or, equivalently: 

Dn{p,q) = M{Dp,Dq). (14) 

4 New approach for the Frechet bounds 

Prom our main result, Theorem l3.lt we obtain Frechet bounds by methods 
which seem to be new to the best knowledge of the author. 

4.1 Upper bound 

The Frechet upper bound is obtained as a direct consequence of Theorem l3.1[ 

Corollary 4.1 (Frechet upper bound) Let p and q be two elements of 
Mat„i(M-|_). Under the condition thatp^^l = q^l the set7i{p,q) is not empty 

and the upper Frechet bound Fmax is such that DFj^axD^ "= -^max is the 
greatest sub-solution of the following max-plus linear system of equations: 

^ (15) 
l^OF =q^D^ ^ ' 

that is: 

F^ax = {{Dp) /I) A {^\{q^D^)). 

Proof. Let us study (jlSp when replacing = by <. Then, by definition 
of / we have: FQl<Dp<^F< {Dp) /I. And by definition of \ we 
have: F < lJ'\{q^D'^). The two previous inequalities are equivalent to: 
F < {{Dp) /I) A {l^\{q'^ D'^)) = Fmax- Now, we have to prove that (A): 
^max Ql = Dp and (B): 1^ F^nax = q^D^. 
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Let us prove (A). 
For all i = 1, 

(Fmax l)i = ®]=i{{Dp)/l) A il^\{q^D^))i,, 1 

TnT\ 



{Dp)i A (e"=i(g^ )j) l^by lattice distributivity 
(Dp), A (q^D^)^ (Vj: (q^D^), < {q^D^^ ) 
(Dp), A (Dp),, ((6^Z)^)„ = g^l = p^l = (Dp) 
{Dp)i ( Vi: {Dp)i < (Dp), 



We prove (B) similarly. Hence the result is now achieved. □ 
4.2 Lower bound 

In this section we obtain Frechet lower bound by a greedy algorithm based 
on max-plus linearity of the Frechet problem and the well-known Monge 
property (see e.g. 0) of elements of the set Dist„„ defined by ([TT]) . that is 
for all F € Dist^n: 

(M). Vi,i = 0, . . . n - 1: FijQ Fi+ij+i > F,j+i F^+ij, 
with the convention that \/k Fq^/^ = Fo,fe = 0- 

Let p,q G Matni(M+) be two vectors such that i = 1 = 1. Let us 

consider the following algorithm. 

Lower(n,p,q) 

Vi = 1, . . . , n, Fj^n ■= {Dp)i ; (a) 
Vj = l,...,n, Fn,j:= {q^D^)j ; (b) 
For j = n — 1 to 1 do 

For z = n — 1 to 1 do 

end 
end. 

Proposition 4.1 The algorithm Lower computes the Lower bound of the 
Frechet contingency array problem. 
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Proof. The initial conditions (a) and (b) of the algorithm Lower come from 
the max-plus linearity of the Frechet problem and that the Monge property 
(M) implies: 

yi<i',j<j', Fi^j<F,,j,. 

The proof is obtained by recurrence (see e.g. the detailed proof of this result 
by Frechet himself [5, p. 13]) 

Denoting ai = {Dp)i, (3k = {q^D'^)k we have to prove that the loop in- 
variant of the algorithm Lower corresponds to the Frechet lower bound, i.e. 

It is easy to see that the previous relation is true for I = n with k = 
1, . . . , n and for / = 1, . . . , n with k = n. Now, let us assume that the loop 
invariant is true for {k,l) > {k,l) / (i, j). We have: 

Fij+i Fi+ij = (1-1 pj+i ©1)0 (1-1 ai+i pj 1) 
= 2-1 Qi ai+i Pj Pj+i 

©1-lQi Pj+l © 1-1 tti+l Pj © 1. 

Because (M+, <) is a totally ordered lattice: 

Fj+ij+i = 1-1 ai+i (3j+i © 1 € {1^1 Ui+i (3j+i, 1} 
Thus, we have two cases to study. 

Irst case : Fj+i j+i = l^i a^+i 
Let us compute: 

Fi+ij+i Fij+i Fi+ij = 1-1 Oi f3j © Oi a-^^ © Pj 

©10 (oi+i 0/3j-+i)-i. 

Then, we just have to remark that: ai a~^\ = p^_^^ < 1, Pj Q Pj'^i = 
qj+i < 1 and ^i+ij+i = 1"^ a^+i pj+i ^ 1 (oi+i Pj+i)'^ < 1. 
Thus, 

= F;^ij+i Fij+i J © 1 = 1-1 ai /3j © 1. 
2nd case : Fj+ij+i = 1 

Fi+ij+i Fi+ij = Fij+i Fi+ij 

= 2-1 Oi+i Pj © l-iftj 
©1-1 Oi+i ,3j © 1 
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which could be rewritten as follows: 



Fij = 1 ^QaiQPjQ{l ^0 a^+i pj+i qj+i pi+i) 1. 

We remark that qj^i,pi+i > 1. We also note that Fj+ij+i = 1 ^ 
aj+i < 1. Thus, we deduce because is non-decreasing and the 
definition of the standard order < that: 

Fij > QaiQ/3j 1. 

On the other hand Fj+ij+i = 1 =^ Pj < qj^i- And we 
deduce that: 

Fij < ^7+1 (1'^ a»+i Pj+i qj+i Pi+i) 1 
= 1"^ QaiQ Pj © qj^-^ © 1 
= 1"^ QaiQPj®! (because g^^^ < 1) 

We conclude because < is antisymmetric. 

□ 



5 Conclusion 

In this paper we proved that the Prechet correlation array problem is max- 
plus linear in the space of cumulative distribution function Dist„n defined 
by (|lip . This remark leads to new methods to obtain the Frechet bounds. 

As a further work we would like to extend results of the paper to the 
continuous case based on the remark that : 

fix, y)dy = sup I / /(x, y)dy 
for all nonnegative functions / such that f{x,y)dy exists. 
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